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SUNMARY 


A method is described from which the interaction of an elastic structure 
with an infinite acoustic fluid can be determined. The displacements of the 
structure and the pressure field of the immediate surrounding fluid are modeled 
by finite elements, and the remaining pressure field of the infinite fluid 
region is given by an analytical expression. This method yields a frequency 
dependent boundary condition for the outer fluid boundary when applied to the 
frequency response of an elastic beam in contact with an acoustic fluid. The 
frequency response of the beam is determined using NASTRAN, and compares 
favorably (1-2% error) to the exact solution which is also presented. The 
effect of the fluid on the response of the structure at low and high frequencies 
is due to added mass and damping characteristics, respectively. 


INTRODUCTION 


The interaction of an acoustic fluid with an elastic solid has received 
considerable attention in the literature. Some areas of investigation in the 
frequency domain include underwater vibrations, vibrations of liquids in elastic 
containers, and the evaluation of the near and far pressure field of an 
acoustical fluid surrounding a sinusoidally excited elastic structure. A finite 
element modeling of the combined problem was formulated by Zienkiewicz and 
Newton (ref. 1) . Their finite element modeling of the displacements of a 
structure and the pressure field of a finite acoustical fluid leads to a 
system of unsymmetric linear equations to be solved. 

Problems involving a finite domain can at least conceptually (and usually 
practically) be modeled using finite elements (see ref. 2, for exanple) , but 
those problems involving an infinite fluid domain must necessarily be modeled 
with only a finite portion of the fluid if the finite element method is to be 
used. The appropriate boundary condition at the truncated fluid boundary is 
often in doubt. Zienkiewicz and Newton (ref. 1) suggest a system of dashpots 
at this outer fluid boundary, but it will be shown that this is the proper 
boundary condition only in the high frequency limit. This paper formulates the 
the boundary condition that should be applied at this outer boundary, and shows 
how this condition is incorporated into the finite element method. To this end, 
the fluid is divided into a region immediately surrounding the structure (which 
is to be modeled by finite elements) and an infinite region. Within the 


infinite region a series expansion is chosen for the pressure, the coefficients 
of which are unknowns of the problem. Hunt et al (ref. 3) have a similar 
model, except the pressure field in the infinite region is given by the sur- 
face Helmholtz integral equation. In any case, the expression for the pressure 
field in the infinite region identically satisfies the governing wave equation 
and the proper boundary conditions at infinity. 

A variational principle, presented specifically for a beam with one face 
in an acoustical fluid, suggests the proper coupling not only between the 
structure and the fluid but also between the finite and infinite fluid regions. 
If a coordinate surface is chosen as the outer fluid boundary, orthogonality 
relationships of the series expansion may be used to satisfy continuity of the 
pressure field at this boundary. This orthogonality allows the coefficients of 
the series expansion to be eliminated as unknowns from the problem, and results 
in an additional stiffness matrix for the nodal pressures on the outer boundary. 
This matrix is full, symmetric, and frequency -dependent, and is inplemented in 
NASTRAN by direct matrix input. 

This method is applied specifically to the frequency response of a sinply 
supported beam with one face in contact with an infinite acoustical fluid (2-D 
problem). The exact solution for the frequency response of the beam is 
presented, and the finite element results compare favorably with the exact 
solution. It is also shown that at low frequencies the effect of the fluid on 
the structure is an added mass, while at high frequencies it is a damping. 
Moreover, the far field pressure in the infinite region can be determined from 
the series expansion once the nodal pressures at the outer fluid boundary are 
known. 

While this method is applied for a 2-D frequency response, it can be 
generalized to the response of a 3-D elastic structure in an infinite acoustic 
fluid. The outer surface of the fluid must be a coordinate surface of a space 
in \diich the wave equation is separable since the orthogonality of the series 
expansion on this surface is used. Once again, the structure and the portion 
of the fluid between the structure and this coordinate surface are modeled by 
finite elements. Unfortunately, the additional stiffness matrix couples all 
the pressure nodes at the outer boundary, and, in general, is frequency 
dependent. If the frequency is specified, the additional stiffness matrix is 
known, although in general it could increase the bandwidth of the problem. On 
the other hand, for determining the submerged natural frequencies of structures 
an iterative procedure is necessary since the natural frequency is unknown. 


A VARIATIONAL PRINCIPLE 


It is convenient in applying the finite element method to have a varia- 
tional principle on which the discretized finite element model can be based. 

Such principles involving the displacements of an elastic structure can be found 
in references 4 and 5; similar principles for fluid mechanics problems are 
presented by Olson in reference 6. Gladwell (refs. 7, 8, 9) presents varia- 
tional theorems for the acoustic fluid for both pressure and displacement 
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formulations. For the coupled structural -fluid problem, a suitable variational 
formulation can be found by properly combining those for an elastic structure 
and an acoustic fluid. Such a principle is a reliable basis and guide for 
numerically solving a fluid- structure problem using finite elements. Moreover, 
with the fluid divided into a finite region (modeled by finite elements) and an 
infinite region (fluid described by an analytical expression) , the variational 
formulation will necessarily point to the proper coupling of each. 


Finite Fluid Region 


A simply supported beam is shown in figure 1 which has one side in contact 
with a finite acoustic fluid and subjected to a sinusoidal load per unit length 
of w(x)ei^t. The deflection of the beam in the y-direction, u(x)ei^b, 
satisfies the differential equation 

A 7 

El — ^ - mft^u = - p(x,0)h + w(x) (1) 

dx^ 


where E is the modulus of elasticity of the beam, I is the moment of inertia, 
m is the mass per unit of length of the beam, and h is the depth of the beam in 
the z-direction. The pressure p(x,y)e^^'t of the fluid region A satisfies the 
wave equation 



A _ 1 
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( 2 ) 


where c is the speed of sound in the fluid. It is also assumed that 


^ = 0 on S (3) 

8n 

where S, shown in figure 1, is the boundary of A excluding the beam's surface. 
On the surface of the beam, it is also necessary to enforce (see refs. 1, 2) 
the condition, which comes from conservation of momentum, that 

= - p ii = p u on y = 0 (4) 

where p is the density of the fluid. 

It is possible to formulate a mixed variational principle that will incor- 
porate both equations (1) and (2) and the appropriate boundary conditions for 
each. Consider the functional F(u,p) given by 
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F(u,p) = y / El( — ^ ) dx - ymQ / u dx + / p(x,0)hudx 
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The first two terms are the strain energy and kinetic energy of the beam, 
respectively. The next two terms are minus the work done by the pressure and 
forcing function on the beam, aitd the last two terms represent the kinetic and 
potential energies of the fluid, respectively. The functional F is a function 
of both the structural displacements u and the fluid pressure p. If independent 
variations of F are taken with respect to u and p, it follows that 

£ 4 

6F(u,p) = / [eI - mP^u + ph - w(x)l 6u dx 

0 L dx^ J 



If u and p are found such that 

6F(u,p) = 0 (7) 

with any trial function u satisfying 

u(0) = u(£) = 0 (8) 

then it can be seen from equation (6) that u and p necessarily satisfy 
equations (1) and (2) and the boundary conditions given by equations fS") and 
(4). 


Coupling of the Infinite Fluid 


If the region of the fluid is infinite, as shown in figure 2, the fluid is 
subdivided: the finite element description of the pressure in the fluid is used 
in a finite region surrounding the structure, and an analytical expression 
(which identically satisfies the wave equation) is used in the remaining 
infinite region A 2 . In order to properly couple the two solutions, the pressure 
field must be continuous and consistent with the variational principle. The 
functional F(u,p) in equation (5) now contains two additional terms which are 
the same as t!ie last two terms but integrated over the remaining infinite 
region. The analytical expression for p in this region satisfies the proper 
boundary conditions at infinity (Somerfeld radiation condition) . Variations of 
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F taken with respect to u and p in both regions give 

5F(u,p) =/ |eI m u + ph - w(x) 6u dx + 6p h ds 

0 L dx -I S 

.El - El 4 U jY^ |£ . ulspdx 

dx Q dx IqOLp^ J 

Aj^ pfi \3x 3y pc y pa s 1 

A-, pa \9x 3y pc / pa s 2 


where s is the boundary between the finite and infinite fluid regions. With 
the analytical function p identically satisfying the wave equation, the 
coefficient of 6p in the next to last term is identically zero. Hence, only 
the term 

( 10 ) 

pa^ s ^^s 

which is the loading of the infinite region on the finite must be included. In 
the integral I, p is given by an analytical expression which must match the 
finite element nodal approximation on the fluid- fluid interface. Assume the 
analytical expression for p is given by an expansion 

N 

p = s 
n=l 

where the Aj^'s are undetermined coefficients and the functions fj^(x,y) 
identically satisfy the wave equation. Equation (11) , together with the finite 
element description of the pressure at the interface and the continuity of the 
pressure field, will permit the integral in equation (10) to be evaluated. The 
continuity of the pressure field can be easily obtained by choosing an outer 
boundary on which the orthogonality of the functions fjj(x,y) can be used. The 
evaluation of equation (11) will be carried out specifically for the frequency 
response of a beam in an infinite fluid. 


FINITE ELEMENT FORMULATION 


by 


Beam and Neighboring Fluid 

The finite element method approximates the displacements u of the beam 
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(12) 


u = E N. u. 
i " ^ 


where u. is a generalized nodal displacement, and Nr is a shape function for 
the displacements o£ the beam. Similarly, the pressure field in the fluid is 
approximated by 


P 


IN^Pi 


(13) 


where pj^ is a nodal pressure. Substituting these approximations into equation 
(6) and interpreting equation (7) to mean that partial derivatives with respect 
to nodal displacements and pressures should equal zero, the following set of 
equations is determined: 


where 
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Miltip lying the second set of equations by (pcf^) gives 
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(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 

(21) 


This form is the same as that derived by Zienkiewicz and Newton (ref. 1) , but 
equations (21) are based on a variational principle. The set of equations (11) 
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can be generated by NASTRAN as outlined in reference 2, although matrices L and 
LT are inputted directly by IMEG cards. 

For more conplicated geometries and structures, the form of equation (21) 
is unchanged. While the variational principle and the finite element formula- 
tion were given specifically for the elastic beam-acoustic fluid problem, they 
can be easily generalized to account for an elastic structure bounded by an 
acoustic fluid. 


Infinite Fluid Coupling Matrix 


The loading of the infinite fluid on the finite portion is found by 
computing the integral 

I --3 

p(2 s 2 

over the outer fluid face (see eq. (10)). This integral is to be discretized 
and then added to the set of equations (14) . At the fluid- fluid interface the 
pressure is given by 

M 

p = E N.p, (23) 

i=l ^ ^ 

where is the shape function for the pressure in the fluid evaluated at the 
fluid-fluid interface, pj^ is a nodal pressure on the face, and M is the number 
of pressure nodes at the fluid- fluid interface. The 6p in equation (22) is 
equal to the partial of p with respect to pi (vdiich is equal to Nj^ from equation 
(23)). Then the term 

i/I^N.hds (24) 

p(2 S 2 

is added to the pj^ equation of equations (14) . 

Consider the frequency response of a beam with one side immersed in an 
infinite acoustic fluid, as shown in figure 3. Both the displacements of the 
beam and the pressure field of the neighboring acoustic fluid are modeled by 
finite elements. The pressure field in the fluid for y>b must be bounded and 
satisfy 

i4.24.-(|)2p C25) 

ax ay^ 

with the bounciary condition that 

p = 0 at X = 0 and x = I (26) 

and the condition that only waves outgoing from the structure are allowed. 
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( 27 ) 


Separation o£ variables leads to the following expression for p: 

N 


with 


_ A rtnrx. ‘“ny 

p= E Asin(-— )e 

n=l ^ ^ 


2 ^ ^m^2 _ n 2 
n ^ a-’ '■c'' 


(28) 


The N arbitrary coefficients are yet to be determined, and ct may be either 
real or imaginary depending on n and Q. “ 

The integral in equation (22) can now be evaluated. From equation (27), 


ay 


9 Ho 


- V A ^nirx-. 

- E A sin (-^)a e 

y=b n=l ^ X. n 


(29) 


Substituting equation (29) into the expression (24) , the following matrix 
expression is added to the left-hand side of equations (14) : 

[G]{A> 




(30) 


where {A} is a vector of the N coefficients A. and [G] is an MxN matrix given 
by J 


-a-jb 




(31) 


The number of unknowns in equation (14) has been increased by N, the number 
of Aj coefficients. An additional set of eqiiations to make the set conplete is 
found by re(quiring the pressure to be continuous; that is, equation (23) must 
match equation (27) evaluated at the interface y=b: 


M N 

'"iPi ■ ^ 


mrx-, ‘“n^ 


i=l ^ ^ n=l 
Multiplying both sides by 


sin(XS) 


and integrating from 0 to with the orthogonality condition that 

0 nj^k 
O' 

J n-k 


( 0 - .,i„/-nTTX., . ,kiTXN j 

J sin(— ) sin(— ) dx = 


(32) 


( 33 ) 


gives 
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{A} = [R]{p} 


(34) 


where [R] is an NxM matrix and is given by 




0 


N. sin (^) dx 


(35) 


The additional N equations from equation (34) form a complete set with equations 
(,14) and (30) . Alternatively, equation (34) can be used to eliminate the 
series coefficients A^ from the expression (30) in favor of the nodal pressures 
of the fluid- fluid interface. Thus, 


[G]{A> ^ [G] [R]{p) 

pq pq 


Define 

[H'l = [G][R] (36) 

with matrices [G] and [R] given by eqmtions (31) and (35) . The effect of the 
infinite fluid on the finite is to add to the fluid stiffness matrix [H] of 
equations (14) the matrix [H' ] . [H' ] is an MxM symmetric matrix which may be 

conplex. It is full, frequency-dependent and couples just the nodes at the 
outer fluid boundary. 


If q-»-0 or equivalently c ->■<», the effect of the infinite fluid on the 
finite fluid is one of stiffiiess. From equation (28), 



a - ->■ ^ as ->■ 0 

(37) 

and 

-a -b 

Gi.-a.e Fy 

(38) 

where 

^ii " / ^i 
13 0 ^ X. 

(39) 

The matrix 

[F] does not depend on the frequency fi. Similarly, 



a^b 

R. . F. . as q->-0 

11 ^ ij 

(40) 

Then matrix 

[H'] = [G] [R] is given by 



“ij ■ °ii ’*''5 

(41) 
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(42) 


Substituting equations (37) , (38) , and (40) into equation (41) gives 




N 

E 

k=l 




vdiere [H* ] is a constant matrix (independent of frequency) and is added 
directly to the stiffiiess matrix [H] of the fluid. If the outer finite element 
boundary is chosen to be the surface of the beam, then the matrices [H] and [Q] 
of equation (14) are zero. Then the second set of equations (14) can be 
written, with the matrix H’ now included, as 


i[H'] {p} + [l'^]{u) = 0 

pq 

This gives 

(p) = -pO^[H’]'^[l'^]{u} 

Substitution of this equation into the first set of equations (14) gives the 
following added mass matrix: 

[M'] = p[L][H’]'^l'^] (43) 

Matrix [M'] is symmetric and full and shows that the effect of the fluid on the 
structure is added mass. 


When or c^O, the effect of the infinite fluid on the finite is a 
pure damping. For, from equation (28) , 


a. 


3 



where i is /TT. Then 




N 

z 

k=l 


''rk ■'sk 


(44) 


[H|] is a pure imaginary matrix, linear in Q, which is to be added to the 
stiffness matrix [H] in equation (14) . Then rewriting equation (14) in the 
form of equation (21) gives 


'K 

.0 


pc 




0 0 "] u 

- 

-0 pcbJ n L- 


u| I 0 0 

+ in 

PI Lo pcBj Ip 


M 

(pc)^l'^ 


0 ■ 

' 

u 

PQ- 

P 


f 

0 


(45) 


where [B] is an MxM damping matrix and an element of matrix [B] is given by 


2 N 

B = 4- E F , F 
rs Z rk rs 
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(46) 


The coefficient defined in equation (39), is equal to (il/2) times the 
Fourier sine coefficient of the shape function Hence the sum of terms on 
the right-hand side of equation (46) is the dot product of the Fourier 
coefficients of the function Nj. with Ng . It can be shown that this dot product 
is equal to (2/£) times the inner product of the shape function N^. with Ng over 
the length of the beam. Thus, 




N^Ngdx 


(47) 


This term is identical to the one suggested by Zienkiewicz and Newton (ref . 1) , 
vdiich is a bomdary condition derived by assuming that the pressure in the 
fluid takes the form of a plane wave. The boundary condition is to be applied 
at a boundary which has been placed "far enough" from the structure and is the 
proper boundary condition only in the high frequency limit. 

If the outer fluid boundary is reduced to that of the beam, then both 
matrices [H] and [Q] in equation (45) are zero. Solving the second set of 
equations (45) for {p} in terms of {u} gives 

{p} = pc(J2i)[B]‘^[L'^]{u) (48) 

Substituting this equation into the first set gives 

[K]{u) + (qi)pcL[B]'^[L^]{u) - q^[M]{u) = {f) (49) 

Then the matrix 

pc[L] 

is a damping matrix, which means that the effect of the fluid on the structure 
at high frequencies (or small c) is danping. 


EXACT SOLUTION - FREQUENQ RESPONSE OF BEAM 


The differential equation of motion for the elastic beam shown in figure 3 
subjected to a uniform load varying sinusoidally in time is 


3x 


dr 


(50) 


vdiere p(x,y,t) is the acoustical pressure which must satisfy the wave equation 



(51) 
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A solution for u of the form 

00 

u(x,t) = E A (52) 

n=l ^ 

is chosen, where the A^^'s are undetermined coefficients of sin , vdiich are 
the out-of-fluid eigenvectors for the singly supported beam. Equation (51) is 
solved using separation of variables. The pressure field is bounded and the 
boundary condition at x=0 and x=i is that p=0. Allowing only outgoing waves 
from the beam leads to the following equation for p: 

00 y ^ 

p - E C sinC-— Je e (53) 

n=l ^ ^ 


where 


2 ,mr>. 2 

ct = (~:r") 
n ^ i ■' 



The Cj^'s are undetermined coefficients to be found by properly coupling the 
fluid and the structure. At the fluid- structure interface, one requires (see, 
for example, ref. 2) 


m = 

9n 



at 


y = 0 


(54) 


Since the fluid and structural modes are uncoupled for this problem, equation 
(54) yields 


c = - H^A 

n n 


(55) 


Substituting equation (55) into equation (53) , and the expressions for u and p 
into equation (50) , gives 

(56) 

with 

2 ^n*iT N 2 2 

“n ' 

When both sides are multiplied by sin (Sj^) and integrated from 0 to t to take 
advantage of orthogonality relationships, the solution is 

u(x,t) = z sin (57) 

n“l , 3 , . . . 
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00 


(58) 


p(x,y,t) = E 

n=l,3. 


“n 


An sin 


nTTX 


”“n^ isit 
e e 


with 




and 


'*'•'0 

* El®'' 
“n 


2 /'TItT'n 2 2 

“n = - y 


(59) 


(60) 


In general may be real or imaginary depending on the driving frequency If 


then for all n, is real. In this case the pressure and displacement are real 
and in phase and no radiation occurs. If > (ttc/ji) , then for some n, Oj. becomes 
imaginary. In this case, both u and p are conplex and out of phase and nence 
radiation may occur. 

Peaks in the frequency response will occur at the in-fluid natural 
frequencies of the beam. Although the in-fluid and in-air mode shape of the 
beam are unchanged in this particular problem, the natural frequency of the beam 
does change. The in- fluid natural frequencies are found by setting the term in 
brackets in equation (56) equal to zero and solving for a. The solution always 
gives < (me/ 1 ) , vdiich means that the in- fluid modal shapes of the beam do not 
radiate . 


RESULTS 


Computations were carried out by NASTRAN using the finite element- 
analytical method previously described. A typical grid is that shown in 
figure 3, where GEAR elements were used to model the beam and 2-D isoparametric 
elements (with quadratic approximation for the pressure) were chosen to model 
the fluid. The usual double numbering of grid points at the fluid- structure 
interface is necessary with this formulation (this procedure is outlined in 
ref. 2) , and the nodal pressures and displacements of the interface are coupled 
through matrix [L] of eqiiation (14) . This matrix is entered into NASTRAN by 
EMIG cards. The frequency dependent matrix [H'] defined in equation (36), which 
models the effect of the infinite fluid on the finite, is also inputted into 
NASTRAN by DMIG cards. The results shown in figures 4 through 8 are for the 
following values: h = 2.54 cm (1 inch), i = 50.8 cm (20 in.), c = 1.460 km/sec 
(5.748x104 in/sec), E = 206.8 GPa (3x10^ psi) , I = 3.468 cm4 (.08333 in4) , 
m = 7.827 g/cm5 (7.324x10-4 Ib-sec2/in4) , p = 1.029 g/cm^ (9.633x10-5 lb-sec2/ 
in4), Wq = 1.751 N/cm (1 lb/ in). The solid line in each of these figures is 
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the exact solution given by equations (57) through (60) . The finite element 
solution is shown at specific plotted points. 

Figure 4 is a plot of the magnitude of the pressure at the center of the 
beam versus the driving frequency u. The pressure peaks at approximately 
1062 rad/sec and 11350 rad/sec, which are the in- fluid natural frequencies of 
the beam for modes n = 1 and n = 3 (these values can be determined by solving 
equation (56) for Q. with wq = 0) . There is a discontinuity in the slope of the 
curve for Q ~ 9029 radians/sec, vdiich is the frequency at which radiation 
occurs (0 = TTc/Jl = 9029 rad/sec). For frequencies greater than uc/l, energy is 
being carried away by the outgoing pressure waves and the beam is said to 
radiate. In this case net work is done by the forcing function. 

Figures 5 and 6 are plots of the magnitude of the beam's displacement at 
its center as a function of the driving frequency Q. In figure 6 the displace- 
ment shows the discontinuity in slope that the pressure exhibits when the beam 
begins to radiate. As JJ goes through ttc/ji, the displacement of the beam 
increases corresponding to the reduction in pressure. 

The variations of the phaise angles of the pressure and displacement with 
frequency are shown in figures 7 and 8, respectively. For Q < -nc/i (9029 
rad/sec), the displacement is in phase (or 180° out of phase) with the driving 
force and no work is done. For > -nc/i, the displacement is out of phase with 
the driving frequency and radiation occurs. The only exception to this 
condition occurs when a approaches a natural frequency. The mode shape for that 
frequency dominates, and, since the in-fluid mode shapes of the beam do not 
radiate, the phase angle of the displacement is in phase (or 180° out of phase) 
with the driving force. 

Figures 4-8 show that the finite element solution obtained through NASTRAN 
was reliable in modeling the elastic beam in the infinite acoustical fluid. 

The errors of the results shown were 1-2% for the grid shown in figure 3. The 
same accuracy was also obtained at a few frequencies in vdiich the outer fluid 
boundary was chosen to be that of the beam (that is , b=0) . In these cases , 
matrix [H'] (defined in eq. (36)) corresponds to the nodal pressures at the 
fluid- structure interface. 

For the limiting case of the effect of the fluid on the structure is 

an added mass; this effect is approximated within the finite element method by 
modeling the structure with NASTRAN and adding to the mass matrix generated by 
NASTRAN the additional mass matrix [M' ] given by equation (43) . The natural 
frequencies and mode shapes of this conputation agreed favorably (less than 1% 
error) with those from the exact solution. The exact solution is determined 
from equation (56) by solving for with wg = 0 and c^«>. 

CONCLUSIONS 

The boundary condition at the truncated fluid boundary of an infinite 
acoustical fluid is, in general, frequency dependent. For a finite element 
formulation this condition leads to a stiffness matrix [H'J which is added to 
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the stiffiiess matrix of the fluid. [H'] is a full, symmetric, complex, 
frequency dependent matrix which couples the infinite region to the finite 
region and involves only the outer boundary nodes. If the driving frequency is 
specified (in the case of frequency response), the coipling matrix [H'] can be 
inputted into NASTRAN by H^IG cards. The conputation of eigenvalues and eigen- 
vectors, on the other hand, would necessarily involve an iteration scheme since 
the frequency of the mode shape is not known. 

Although only the portion of the fluid immediately surrounding the 
structure is modeled by finite elements, the infinite fluid region is 
effectively modeled through the coupling matrix [H*]. Moreover, the far field 
pressure can be determined once the outer boundary pressures are conputed. 

This pressure is given by eqxoation (27) with the series coefficients {A} 
determined from equation (34) . If a finite portion of the fluid is modeled 
without including the boundary condition matrix [H'], then the fluid region is 
actually a finite domain. Not only is it inpossible to determine the far field 
pressure but also some of the eigenvalues and eigenvectors found can be shown 
to be associated with the finite problem. These additional modal values do not 
appear if [H*] is included. 

This type of finite element-analytical solution, presented here for a 
two-dimensional problem, can be readily generalized to a three-dimensional 
problem of modeling an elastic structure in an infinite acoustic medium. In 
this case the outer fluid boundary would be a sphere, and the pressure field 
would be given by an expansion of spherical harmonics. The frequency dependent 
matrix could be generated in NASTRAN by program modifications and would be 
accessed through IMVP alters. IMfortunately, because of the full coipling of 
all the outer boundary nodes, the increase in the bandwidth for a three- 
dimensional problem might make the conputer cost prohibitive. 
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INFINITE FLUID REGION 


Fluid Element 


Figure 3. Finite Element Grid Showing Both the Fluid and Beam 
Elements. Nodal Pressures of the Fluid and Nodal Dis- 
placements of the Structure are the Unknowns. 



Figure 4. Magnitude of the Acoustical Pressure at Point A (Center of Beam) 
as a Function of the Driving Frequency, Peaks Occur at the 
In-Fluid Natural Frequencies Given by Equation (56) . 
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Figure 6. Magnitude o£ the Normal Displacement of the Center Beam Point (Pt. A) 
a Function of the Driving Frequency. Peaks Occur at the Natural Frequencies 

Given iDy Equation (56) . 
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NASTRAN Solution 
— Exact Solution [Eq. (58)] 


DRIVING FREQUENCY (radi^ns/sec) 


Figure 7. Phase Angle o£ the Acoustical Pressure at the Center o£ 
the Beam as a Function o£ the Driving Frequency 



X NASTRAN Solution 

— Exact Solution 
[Eq. (57)] 



DRIVING FREQUENCY (radians/sec) 

Figure 8. Phase Angle o£ the Normal Displacement o£ the Center 
Beam Point as a Function o£ the Driving Frequency Q. 



